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The remarkable properties of the symbiotic star AE Circinus 
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ABSTRACT 

We present new optical spectroscopy and photometry, 2MASS infrared observations and 
24 years of combined AAVSO and AFOEV photometry of the symbiotic star candidate se. 
The long-term light curve is characterized by outbursts lasting several years and having a slow 
decline of 2 x lO^** mag/day. The whole range of variability of the star in the V band is 
about 4 magnitudes. The periodogram of the photometric data reveals strong signals at ~ 342 
and 171 days. The presence of the emission feature at A 6830 A at minimum and the detection 
of absorption lines of a ~ K5 type star confirm the symbiotic classification and suggest that 
AECir is a new member of the small group of s-type yellow symbiotic stars. We estimate 
a distance of 9.4 kpc. Our spectrum taken at the high state shows a much flatter spectral 
energy distribution, the disappearance of the A 6830 A emission feature and the weakness of 
the Hell 4686 emission relative to the Balmer emission lines. Our observations indicate the 
presence of emission line flickering in time scales of minutes in 2001. The peculiar character 
of ae is revealed in the visibility of the secondary star at the high and low state, the light curve 
resembling a dwarf nova superoutburst and the relatively short low states. The data are hard 
to reconciliate with standard models for symbiotic star outbursts. 
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1 INTRODUCTION 

AECir (S32 = StHA32, Q2000 = 14:44:52.0, <520oo = -69:23:35.9) 
was classified as a RGB star in the Fourth General Catalog of Vari- 
able Stars. This classification was rejected by Kilkenny (1989), 
who found strong emission lines of H I and weaker emission lines 
of He II and other species in one spectrogram spanning a wave- 
length range of 3400-5 100 A. Kilkenny noted the lack of forbidden 
lines and the anomalous strong He II 4686 emission for this object, 
and suggested a symbiotic classification, although no lines of the 
cool component were detected. He also mentions the photometric 
variability of the star, ranging from 12 to 14 mag. (from visual esti- 
mates) in hundreds of days. B, V photoelectric measures by Law- 
son & Gottrell (1990) taken in an interval of 107 days show the star 
with V magnitude between 13.52 and 14.51 and color B — V be- 
tween 0.95 and 1.43, the object being redder when fainter. Based 
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on the spectroscopy by Kilkenny, AE Cir was listed as a suspected 
symbiotic in the catalog of Belczynski et al. (2000). Symbiotic stars 
have been reviewed recently by Mikolajewska (2007). In this paper 
we present new photometry and low and high-resolution spectra of 
AE Cir and investigate their emission line properties and symbiotic 
nature. 

In Section 2 we give details of our spectroscopic observa- 
tions and show the methods used in data reduction and analysis. 
In that section we also analyze available visual long-term photo- 
metric records along with our own CCD photometry. In Section 3 
we analyze our data, giving the main results of our research. In Sec- 
tion 4 we discuss our results in the context of symbiotic stars and 
present our conclusions in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 

In this section we present photometric time series obtained at La 
Silla and Cerro Tololo Inter American Observatory (CTIO) dur- 
ing several seasons. We also present medium resolution spectra of 
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ae obtained in May 29 and 30, 2000 with the 1.54 m Danish tele- 
scope at La Silla European Southern Observatory (ESO), high res- 
olution spectra obtained at the Paranal Observatory with the ESO- 
VLT-UT2 telescope during March 16, 2001 and a medium resolu- 
tion spectrum obtained at the CTIO with the 4m Blanco telescope 
during May 04, 2007. Instrumental setups and observational details 
are given for all observing runs in this section. The spectroscopic 
observing log is presented in Table 1. 

2.1 ESO and CTIO photometry 

AECir was observed with the Dutch 0.9m telescope at La Silla 
ESO observatory in the night of 5/6 April 2001 in the V band to 
check for short-term variability. In a total of 82 images with 150 
sec exposure time each, AECir is seen at constant brightness of 
V = 13.65 ± O.IO with a scatter of ± 0.03 mag, consistent with 
statistical fluctuations. 

We also observed AECir with the CTIO 1. 3m telescope op- 
erated under the SMARTS consortium in the seasons 2004 and 
2005. Observations were done in pairs of two exposures in B and R 
each, on a total of 39 nights. Relative photometry was obtained us- 
ing daophot in IRAlQ Absolute photometric calibration was done 
against the USNO-Bl star 0205-0393810 with B = 15.0 and R = 
13.9 mag, respectively. Our relative photometry is accurate to 0.05 
mag, as verified against a number of nearby constant stars. 



(about 0.06 A rms, or 4 km s^^ at H/3). The atmospheric absorp- 
tion bands in the red spectra were removed from the spectra us- 
ing the "telluric" IRAF task. For that we used a template obtained 
by normalizing the standard star spectrum to the continuum with 
a high order polynomial. This high-order function was successful 
in fitting not only the continuum, but also the relatively broad stel- 
lar lines of the template without affecting the narrow telluric lines. 
This method allowed a good removal of telluric lines in our red 
spectra. A flux calibration was also performed with the pipeline, 
with the zero point accurate at the level of 10%. A correction factor 
to the photon counts was applied to the spectra to account for the fi- 
nite slit width and the seeing disc size reported in the ESO Ambient 
Conditions Database. 

2.4 CTIO medium resolution spectroscopy 

We obtained one spectrum of AE Cir on May 04, 2007 with the RC 
spectrograph mounted in the CTIO 4m Blanco telescope at a reso- 
lution of 4.82 A. Reductions were done in the usual manner with 
IRAF, calibrating in wavelength with HeNeAr lamps that provided 
62 lines that were fit with a 6th order Chebychev function resulting 
in a rms of 0.24 A (1 1 km s^^ at Ha). We also flux calibrated our 
spectrum with the aid of the standard star LTT 4816 observed at the 
beginning of the night. No telluric correction was applied for this 
spectrum. 



2.2 ESO DFOSC low resolution spectroscopy 

The Danish Faint Object Spectrograph and Camera (DFOSC) was 
used with several grisms yielding a combined wavelength range 
of AA 3500-9500 A. A slit width of 1.5"was chosen in order to 
match the typical point spread function at the focal plane of the 
telescope. This resulted in spectral resolutions of 4 A (grism 7), 
7.5 A (grism 5) and 14 A (grism 15). The instrument was rotated 
before exposures to ensure optimal orientation relative to the paral- 
lactic angle. Reductions were done in the usual manner with IRAF. 
The spectra were flux-calibrated with the standard stars LTT 3864 
and LTT 7987 observed with slit width of 5"and were calibrated 
in wavelength with He-Ne exposures taken typically after half-an- 
hour sequences. A typical rms of the wavelength calibration func- 
tion for grism 7 was 0.1 A (5 km s~^ at Ha) and the zero point for 
the flux calibration, corrected for the finite slit width, is accurate at 
the level of 15%. 



2.3 VLT UVES high resolution spectroscopy 

We used the high-resolution optical spectrograph of the VLT, 
UVES, in visitor mode. Our setup and the slit of 0.9"yielded a re- 
solving power of ~ 50.000. We obtained spectra in the blue spec- 
tral range using the blue arm and cross disperser CD#2 with cen- 
tral wavelength A 4630 A, and red spectra using the red mode and 
cross disperser CD#3 centered at A 6100 A. The two CCDs in the 
red camera provided two quasi-simultaneous spectra separated by 
a gap centered on the central wavelength. Reductions were made 
using the UVES pipeline software, providing wavelength calibra- 
tions - based on fits of Th-Ar lines - with an accuracy of 1.6 pixels 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the NSF. 



3 RESULTS 

In this section we present and discuss photometry obtained by 
the American Association of Variable Star Observers (AAVSO) 
and the Association Francaise des Observateurs d'Etoiles Vari- 
ables (AFOEV) that is available at |http://www.aavso.org/l and 
http://cdsweb.u-strasbg.fr/afoev/, respectively. Our own photome- 
try is also discussed along these data. We also give a general de- 
scription of the spectra of our target and discuss the published 
infrared photometry. We show evidence for a cool stellar compo- 
nent and determine its spectral type. We also provide an estimate 
of the reddening and distance and present spectroscopic measure- 
ments for the emission lines. The emission line ratios are compared 
with model predictions and physical conditions of the line emitting 
region are estimated. 

3.1 AAVSO, AFOEV and CTIO long-term photometry 

24 years of visual photometry are summarized in Fig. 1 . In spite of 
the limited accuracy of visual estimates (formally ± 0. 1 mag, but 
probably ± 0.2 or 0.3 mag), the general long-term pattern of vari- 
ability appears clear. After a prolonged epoch of slow decline from 
maximum, the star rapidly enters in a faint stage, reaching magni- 
tudes usually ~ 3 mag fainter than those observed at maximum. 
The minimum lasts ~ 900 days and is characterized by irregular 
short-term variability. Then the star recovers its maximum bright- 
ness and the cycle starts again. The prominent maxima observed 
at HJD 2 447 759 and 2 452 390 indicate a time scale of recurrence 
of about 4630 days and the whole range of variability of the star 
in the V band is about 4 magnitudes. However, this pattern does 
not repeat strictly periodically; the duration of the minimum and 
the subtle characteristics of the light curve are different within the 
few cycles covered by the data. We note that the star shows larger 
amplitude variability during the broad minimum than in maximum. 
We also note that the main outburst covered by the data could have 
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Table 1. Log of spectroscopic obsei-vations. N is tlie number of spectra. We give the resolving power R, HJD (-2450000) at the start and 
the end of the observations and the mean seeing value. 



Date(UT) 


Observatory 


N 


HJD{start) 


HJD(end) 


Grism/mode 


Range(A) 


R 


exptime (s) 


slit/seeing (") 


29/05/00 


La Silla 


1 


1694.7079 


1694.7079 


15 


3500-9500 


400 


300 


1.5/1.5 


29/05/00 


La Silla 


24 


1694.7152 


1694.7724 


7 


3840-6820 


1500 


120 


1.5/1.5 


30/05/00 


La Silla 


64 


1695.6748 


1695.8160 


7 


3840-6820 


1500 


120 


1.5/1.5 


30/05/00 


La Silla 


5 


1695.8226 


1695.8305 


5 


5080-9370 


800 


120 


0.9/1.8 


16/03/01 


Paranal 


28 


1984.8623 


1984.9103 


red 


5175-6075 


50000 


120 


0.9/1.8 


16/03/01 


Paranal 


28 


1984.8623 


1984.9103 


red 


6119-7105 


50000 


120 


0.9/1.8 


16/03/01 


Paranal 


15 


1984.8224 


1984.8519 


blue 


4014-5241 


50000 


150 


0.9/1.8 


04/04/07 


Tololo 


1 


4194.5472 


4194.5472 


KPGL3-1 


3700-7230 


2000 


900 


1.0/NA 
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Figure 1. Above: The AAVSO (asterisks) and AFOEV (dots) light curve of 
AECir along with the Lawson & Cottrell (1990) magnitudes (stars), our V 
magnitude obtained at La Silla in 2001 (filled square) and our spectropho- 
tometric V magnitude of year 2000 (filled circle). We indicate times for 
2MASS, DENIS, ROSAT and our 2001 and 2007 spectroscopic observa- 
tions. Below: Smoothed AFOEV data after subtracting the long-term ten- 
dency showing variability in time scales of months. 



been started rapidly, in just 16 days. This rapid rising time is based 
on two datapoints, at HJD 2447 614 and HJD 2447 684. Another 
observation at HJD 2447 708, i.e. already inside the outburst, also 
indicate a magnitude around 14.8 and if true should indicate a short 
fading episode. The outburst rising time will be relevant for the 
interpretation in Section 4. On the other hand the jump to the out- 
burst with maximum around 2 452 390 seems to be preceded by a 
gradual increase of brightness. This could be the rule for AECir 
if we neglect the three mentioned datapoints around the beginning 
of the main outburst. On the other hand, we doubt about the real- 
ity of the maxima at HJD 2 446 892, 2 446 905 and 2 446 9 1 3 , since 
they indicate V about 2 mag. brighter than the average measured 
from nearby observations. Both outbursts show a slow decline of 
~ 2 X 10~* mag/day. We observe also the final part of a third 
outburst at the beginning of the time series. Our V spectrophoto- 
metric magnitude of year 2000, the ESO V average magnitude of 
year 2001, and the times for our spectroscopic observations are also 
shown in Fig. 1, indicating that we observed the star at minimum in 
years 2000 and 2001 and during the high state in year 2007. 

We analyzed the 5-point average smoothed AAVSO light 
curve during the first long decline between HJD 2447 700 and 
HJD 2450 900. The long-term decline was subtracted removing 



the best third order polynomial fit from the data. We analyzed 
residuals with the pdm IRAF package (Stellingwerf 1978) finding 
a period of 342 ± 15 days. A second periodicity of 171 days is also 
observed in the periodogram but with less power. The Period04 
program first developed by Martin Sperl (Lenz & Breger 2005, 
|http://www.astro.uiuvie.ac.at/~dsn/dsn/Period04/l, based on the 
Fast Fourier Transform, yields similar results, although it favors the 
shorter periodicity, probably since it is more sensible to sinusoidal 
variability than programs based on phase dispersion minimization 
algorithms. We argue later that the 342 day periodicity provides 
a more realistic light curve. We note that the periods are highly 
significant; the Fourier periodogram shows the short and the long 
periodicities at the level of lOcr and 8cr, respectively, where a is 
the noise amplitude (Fig. 2). A t-test applied to 10 consecutive 
data bins in the phase diagram shows inconsistency with a random 
distribution at the 0.05 level for all except one case. A visual 
inspection of the data during minima barely reveal the presence 
of any periodicity. Some prominent 1 -point minima observed in 
the light curve do not follow any periodicity (e.g. those at JD 
2453 096, 2453 436 and 2453 835). When analyzing the data 
using 3 consecutive sets of 1200 days length each, we obtain 
periods of 328, 364 and 364 days, respectively. Hence some 
period increase is possible during decline, but the uncertainties 
(several tens of days) are too large for this increase to be measured. 
The few times of maxima detected before the maximum around 
HJD 2 447 700 are coherent with those of maxima after that 
maximum, and the amplitude decreases as the outburst evolves. 
In principle the period close to 1 year could reveal a seasonal 
fiducial periodicity, but the large amplitude observed in the light 
curve lets this appear unlikely. The light curve phased with the 
342-day period is very noisy (Fig. 2), probably a result of the 
large uncertainties of the visual estimates. However, we can see 
an ellipsoidal-like variability characterized by two minima of 
different depth. It is in principle possible that the 342 day period is 
the orbital period of a binary star. In this case, the ephemeris for 
the main minimum is: 

T(HJD) = 2 447 327.3(±3.4) + 342(±15) x N (1) 

We note that the symbiotic nova VI 329 Cyg shows a similar os- 
cillatory behavior, with a period of about 955 days (Chochol & 
Wilson 2001). In that system, polarimetry and radial velocity mea- 
surements prove this period to be the orbital period of the binary 
system. On the other had, the symbiotic star Z And also shows 
a similar behavior during decline of outburst (Formiggini & Lei- 
bowitz 1994), but in this case the oscillations are quasi-periodic 
and shorter than the orbital period. 

The CTIO long-term photometric data are shown in Fig. 3 
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Figure 2. Upper graph: Fourier spectrogram for the AFOEV data taken dur- 
ing tlie long decline, after smoothing and removing the long-term tendency. 
The main peaks are around 342 and 171 days. Lower graph: The same data 
phased with the period 342 days and the 100-point average curve. Stars 
represent the average-subtracted Lawson & Cottrell (1990) data. 



(filled symbols), with the relevant AFOEV data overplotted. We 
clearly resolve the short-term cycles, and find that the color B — R 
changes systematically during the phases of minimum light. This 
variability is not resolved in the AFOEV data, indicating that the 
uncertainty of the visual estimates is indeed larger than the formal 
value of ± 0.1 mag. We observe an eclipse-like episode around 
HJD 2453 620 and an additional smoother - but not completely re- 
solved - minimum about 340 days before, at HJD 2 453 280. This is 
about the periodicity found in the visual data. We cannot compare 
these timings with the predictions of Eq. 1, since the uncertainties 
involved in the period and zero point are too large. Note that our 
mean B — R = 1.0 mag, as compared to that derived from USNO- 
A2 of i? — i? = 0.1 mag which is from non-simultaneous obser- 
vations. We discuss the characteristics of the light curve in Section 
4. 
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Figure 3. Light curve obtained with the 1.3m CTIO telescope (filled sym- 
bols) along with shifted AFOEV data (triangles). Note the eclipse-like 
episode around HJD 2 453 620. 



The magnitudes calculated from the spectra after convolution 
with a Johnson-V filter transmission curve are V = 14.6 in year 
2000, V= 15.9 in 2001 and V= 13.04 in 2007. These magnitudes 
are roughly consistent with the general pattern of the visual light 
curve indicating that the star was observed at minimum during 
2000-2001 and declining from maximum at 2007. We note that our 
V magnitude of 2001 is quite low, the faintest ever observed. The 
average V magnitude 13 and 20 days after our 2001 spectroscopy 
are V= 13.25 and V= 13.65, respectively, indicating a possible 
jump of more than 2 magnitudes in a few days. However, due to 
the large seeing during the 2001 observations, is possible that the 
slit loss was larger than our correction, and therefore we consider 
hereafter our V magnitude value as a lower limit only. Regarding 
the 2007 magnitude, it fits better the tendency of the light curve in 
Fig. 1, but we note that although the star was observed at the be- 
ginning of the night, when the sky was clear, clouds were reported 
1.5 hours after finishing the AE Cir observation, hence we give low 
confidence to our 2007 spectrophotometric magnitude. 



3.2 The average spectrum and spectrophotometric 
magnitudes 

The averaged spectra are shown in Fig. 4. We observe narrow emis- 
sion lines of the Balmer series and also He I and He II emissions. 
The emission bump at 6830 A typical of symbiotic stars is visible 
at minimum in 2000 and 2001, but not at maximum in 2007. No 
forbidden lines are observed. We note high excitation lines as the 
very strong Hell 4868 emission and the presence of the CIII-NIII 
complex at 4630-4650 A. In 2000 and 2001 the stronger line in the 
complex is N III 4640 A, but in 2007 the set of C III and O II lines 
around 4648 Ais the dominant part. The continuum slowly rises to 
the red in 2000 and 2001 but is quite flat in 2007. The NaD absorp- 
tion line is strong and visible in all epochs. We note TiO bands of 
a cool star; these will be analyzed in the next section. We also note 
that during 2007 the Oil 4414 A line is stronger than in previous 
epochs. In 2001 we observe a larger emission strength of He 4686 
compared to year 2000, and a reversal of the 4686/Hq ratio. In 2007 
the object is very bright, and this coincides with a smaller emission 
strength for all lines, especially He II 4686, that is weaker than H/3. 



3.3 Detection of the cool stellar component 

Absorption lines of the cool star are clearly visible in our spectra 
(Fig. 5). We used the indices defined by Kenyon and Fernandez- 
Castro (1987) to find the spectral type of the cool stellar compo- 
nent. 

We measured these indices in the 2000 average spectra for ev- 
ery grating, and obtained [TiO]i= 0.041 ± 0.003 and [TiO]2= 
-0.004 ± 0.009. Then, using the formula given by Kenyon & 
Fernandez-Castro (1987) we found a spectral type ~ K5 consis- 
tent with both indices. The same result is obtained with the 2007 
spectrum, where we found [TiO]i= 0.053. 

In order to check if the cool star shows the same spectral type 
in the spectra obtained in 2001, we subtracted a grid of normalized 
UVES G8— M6 high resolution stellar templates with determined 
MK spectral types (http://www.sc.eso.org/santiago/uvespop/) from 
a region basically depleted of emission lines. The templates were 
allowed to vary in intensity by a scale factor and in wavelength 
by a wavelength shift. The process was performed with the IRAF 
task "telluric", yielding the lowest rms at the best stellar template 




Figure 4. The average spectrum of AECir during the observations of year 2007 (upper frames), 2001 (middle frames) and 2000 (lower frames). The spectra 
are separated in two sections and cut in the y-axis in order to increase the visualization of the spectral features. Some emission lines and the TiO band around 
A 6300 A are labeled for illustration. In the 2007 spectrum, the symbol "x" indicates bad column regions. 
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compatible with tiie cool star spectrum. Our grid included giants 
with spectral types G8, KO, K3, K5, K7, MO, Ml, M4 and M6. 
The best solution was obtained with the K7 III template HD 89736. 
However, the templates in the range KB to Ml also gave a reason- 
able fit, differing with the K7 solution only marginally. The method 
was applied in the regions 5845-6000 A, and 6380-6460 A, obtain- 
ing consistent results in both regions (Fig. 6). The luminosity of the 
cool star cannot be established from our spectra due to contamina- 
tion by line emission of the spectral indices usually employed as 
luminosity indicators in late type stars (e.g. Malyuto, Oestreicher 
& Schmidt-Kaler 1997, Schmidt-Kaler & Oestreicher 1998). 

In order to determine the velocity of the cool stel- 
lar component, we compared our high resolution spec- 
trum in the region 5900-5970 A with a synthetic spec- 
trum of similar resolution calculated with the spectral 
synthesis program SPECTRUM by Richard O. Gray 
( http://www.phys.appstate.edu/spectrum/spectrum.html , Gray 
& Corbally 1994). We run SPECTRUM with a Castelli-Kurucz at- 
mosphere model ( http://wwwuser.oat.ts.astro.it/castelli/grids.html I 
with Teff = 3750 K, logg = 3.0 and solar composition. By 
shifting the object spectrum in velocity until matching the spectral 
lines of the synthetic spectrum we calculated a velocity, with 
respect to the Local Standard of Rest, of Vsec= -25 ± 3 km s^^ in 
2001. For the spectra of years 2000 and 2007 the resolution is too 
low to obtain reliable results with this method. For these spectra 
we used the lines indicated in Fig. 5 to calculate Vsec= -4 ± 7 km 
s"^ in 2000 and Vsec= +2 ± 25 km s"^ in 2007. 

We observe that the NaD(l-l-2) equivalent width is variable: 
2.50, 3.72 and 1.80 A in 2000, 2001 and 2007, respectively The 
values are lower than those expected for a K5 star, ~ 5 A (Tripic- 
chio et al. 1997). This could mean that we never observe the cool gi- 
ant completely, but veiled by some additional light in optical wave- 
lengths. This is more evident in the high state, where the flux is 
flatter and the NaD line weaker. 

We investigated the infrared photometry of £e extracted from 
the 2MASS catalog using the VizieR on-line query form (Ochsen- 
bein et al. 2000). The 2MASS photometry was obtained at JD 
2451 624.7030, i.e. at epoch of minimum light (Fig. 1). It gives J 
= 10.88 ± 0.02, H = 10.05 ± 0.02 and K = 9.74 ± 0.02, hence the 
object is much brighter in infrared wavelengths than in the optical 
region at minimum. The DENIS catalog gives I = 12.42 ± 0.02, J 
= 1 1.01 ± 0.06 and K = 9.82 ±0.05 obtained at JD 2451 274.7594, 
in general agreement with the 2MASS data. DENIS also provides 
the USNOA2.0 magnitudes of 7? = 12.7 and B = 12.8 obtained at 
different epoch. The 2MASS infrared colors J ~ H = 0.83 ± 0.03 
and H — K = 0.31 ± 0.03 are typical for S-type symbiotics and 
probably exclude the presence of dust emission at this epoch (Allen 
1982). The colors of S-type symbiotics have been interpreted as due 
to the presence of a late star, typically an M giant, as opposed to 
D-type symbiotics which contain heavily dust-reddened Mira vari- 
ables. 

We do not have a determination of the color excess at the 
epoch of the infrared observations. However, assuming that AE Cir 
is beyond the galactic dust disc, it should has at least E{B — V) 
= 0.26, the value obtained from the galactic extinction map by 
Schlegel, Finkbeiner & Davis (1998) in the direction of AE Cir (In 
= 312.67, bii = -8.69). This means dereddened colors of (J — H)o 
< 0.76 ± 0.03 and {H - K)o < 0.27 ± 0.03. The first color is 
consistent at the low limit with the value expected for a K5 giant 
({J —H) = 0.73) but the second color is a bit redder compared with 
the expected (H — K) =0.18. The above reference colors are from 
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Figure S. The normalized spectra showing secondary star features. The TiO 
absorption bands used in the determination of the spectral type are labeled. 
Arrows indicate lines used for measuring the cool star radial velocity. Spec- 
tra of two giants of known spectral type taken from the UVES database (see 
text) are shown for comparison. No telluric correction has been applied for 
the science spectra, and UVES spectra were resolution degraded for bet- 
ter comparison. Bad column regions of the 2007 spectrum were continuum 
interpolated. 
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Figure 6. Illustration of the process of searching for the spectral type of the 
secondary star in spectra of March 2001. We show the spectrum of AECir 
in a region depleted of emission lines, along with the spectra of Ml III and 
K7 III template stars, and the result of the optimal subtraction of absorption 
features. 



Koornneef (1983), after transforming to the 2MASS system using 
the equations by Carpenter (2001). 

An estimate for the distance can be obtained assuming maxi- 
mum visibility of the red giant at the minimum ofV ~ 15.5 and 
using reddening E{B — V) = 0.26. We obtain d = 9.4 kpc assuming 
Mv = -0.2 for a K5 giant. 
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Table 2. The emission lines observed in the spectra of ae. We give the equivalent width EW, the full width at half maximum GFWHM, the radial velocity at the line 
center with respect to the Local Standard of Rest, the ionization potential Xion and the maximum line intensity normalized to the continuum. The equivalent width of 
the NaDl and NaD2 lines is the total equivalent width of the three components. All parameters refer to 2001 spectra, unless specified, "nr" means line not resolved, "a" 
asymmetrical line and "b" blended line. 
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Table 3. Table 2 - Continuation. 
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3.4 The high resolution spectrum of year 2001 

We searched for emission lines in the averaged spectra and fit them 
with single Gaussian functions or, in the case of blended fines, witfi 
multiple Gaussian functions with variable central wavelength Ac 
and full width at half maximum GFWHM. Some lines were not 
separated at our instrumental resolution, for instance C II 4267.02 
and C II 4267.27. In these cases a single gaussian fit was performed 
and then the average rest wavelength used as the central wavelength 
for the radial velocity calculation. The emission bump at A 6830 A, 
usually observed in symbiotic stars, is difficult to identify with a 
single line, their GFWHM in 2001 was 26 A, and its EW was 8 
A in 2000, 5.5 A in 2001 and zero in 2007. We find no trace of the 
A 7082 A bump sometimes observed in symbiotic stars. From the 
parameter Ac we measured the radial velocity for every line with 
respect to the Local Standard of Rest (LSR). Results of our analy- 
sis are shown in Table 2. We note that besides HI, Hel and Hell 
emission lines, we find emission of Mg II, Fe II, O II, N III, C II and 
C m. The C II lines detected at A 6576.45 A and A 6581.25 A are 
doubles and show a half peak separation of 33 ± 5 and 26 ± 5 km 
s~^ respectively. The Hues with maximimi ionization potential ob- 
served in the spectra correspond to those of He 11 (54,4 eV). For the 
2007 spectrum, the average radial velocity for 8 emission lines is 
vlsr= -44 ± 8 km s~\ 

We show in Fig. 7 the main averaged emission lines during 
2001 along with the standard deviation obtained after the averaging 
process. We have averaged continuum normalized spectra. In Fig. 7 
we note the asymmetry of He II 4686 and Balmer lines - their max- 
ima are displaced to the red. The Balmer lines are very interesting. 
We note that the asymmetry could be due to blending with Hell 
lines. However, equivalent widths for isolated He II lines are of the 
order of 2 or 3 A whereas those obtained for He 11 lines, if they 
are assumed to be blended with symmetric Balmer lines, are larger 
(Table 2). We also note that in spite of the changes in emission line 
strength, the ratios between the red and violet peak of the Balmer 
lines are very smiilar (2.7, 2.8, 2.6 and 2.3 for H5, H7, H/? and Ha 
lines, respectively) and that this tendency caimot be explained by 



blending with He II lines. Here we have interpreted the shoulders 
observed in the violet part of the Balmer lines as the "violet peak" 
whereas the main peak is the "red peak". We conclude that the 
Balmer lines are intrinsically asymmetric and probably have a very 
small blending with Hell lines. Therefore the parameters for the 
He II lines located near the Balmer lines should be taken with cau- 
tion. Consequently, only isolated He II lines are considered in our 
subsequent analysis. When using the maximum intensity as the ref- 
erence for the radial velocity of the Balmer lines, we obtain values 
that are -10, -19, -20, and km s^^ different from the values listed 
in Table 2, for US, H7, H/3 and Ha. On the other hand, when us- 
ing the line barycenter, we obtained values differing by +3, -i-3, -i-12 
and -1-43 km s"^ from the values listed in Table 2 for the same lines. 
This exercise illustrates the influence of the Balmer line asymmetry 
on the radial velocity measurements. We note that the peak sepa- 
ration is about 150 km for all Balmer hues. In addition, the 
Ha emission wings extends beyond 500 km s~^ from the velocity 
associated with maximum intensity, and shows several structures. 
The red wing has a steep slope and reveals a small but significant 
absorption feature at A 6564.24 A. On the other hand, the blue wing 
has a much less steep slope, and shows a broad bump redwards A 
6560 A. The FWHM measured with the cursor is 2.86 A (130 
km s~^). This is significantly lower than the values obtained with 
the averaged grating-7 spectra of 2000, namely 365 km s~^ and the 
value obtained in 2007, viz. 552 km s~^. Even considering instru- 
mental broadening the Ha line is broader in year 2007. Contrary 
to other epochs, the 2007 Ha emission profile looks symmetrical 
and double peaked with full peak separation of 332 km s^^. We 
note that when compared with our spectra, the (non flux calibrated) 
spectrum shown by Kilkenny (1989) reveals larger emission: this 
author reports equivalent widths of -81 and -55 A for the emission 
lines He II 4686 and H/3, respectively. 

The rms curves in Fig. 7 reflect the net line variability if the 
continuum remained constant during our short observations. Vari- 
ability is not observed in the NaD lines, neither in the bump at 
A 6830 A. However, we do observe variabiUty in others fines. In 
general, this variability is larger in the line center. Maximum vari- 
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Figure 7. Average normalized emission lines (upper profiles, labels at the left axis) and standard deviation (lower profiles, labels at the right axis) during 
March 200 1 . Lines of the C III-N III complex are labeled for illustration. Note that variability is higher in some lines, and mainly restricted to the line center. 



ability is observed in the He II 4686 line, amounting to 30% of the 3.5 The A 6830 A emission bump 

hne intensity. H/3 is variable at the 15% level and Ha at the 2% 

level. 



During the 2001 observing session we do not observe radial 
velocity changes in any of the lines. The He II 4686 line, for in- 
stance, remains at the same position during 45 minutes with a mean 
of 4684.697 ± 0.003 A (rms, 15 spectra). For this reason, we con- 
centrated our analysis on the average spectrum, and measured the 
radial velocity of their lines. The radial velocity distribution of the 
lines have a broad maximum between -63 and -23 km s^^. The 
average is -42 ± 19 km s~^ (rms, compare with -26 ± 20 km 
s"^ in 2000 and -44 ± 8 km s"^ in 2007). The radial velocity 
shows a slight decrease at lines of high ionization potential. A lin- 
ear least squares fit gives v (km s~^) = -30.9(±6.7) - 0.59(±0.22) 
X, with a correlation coefficient of -0.68, where x is given in eV. 
The GFWHM shows the contrary effect; it increases with the ion- 
ization potential. A linear least squares fit, excluding the H I outlier, 
gives GFWHM (km s"^) = H.1{±1.1) + 0.64(±0.24) x, with a 
correlation coefficient of 0.69. In this case the H I lines are clearly 
broader and deviate from the general tendency, a behavior already 
observed in Fig. 7. Since we observe this tendency in all Balmer 
lines, it cannot be due to the phenomenon of Raman scattering, that 
has been suggested to explain the large Ha wings observed in some 
symbiotic stars (Nussbaumer, Schmid & Vogel 1989). 



The presence of the A 6830 A emission feature in AECir is a 
strong argument in favor of the symbiotic classification. At present, 
this feature has not been observed in any other kind of astrophys- 
ical object (Allen 1980). Schmid (1989) proposed that the A 6830 
A emission observed in symbiotic stellar systems is the product 
of resonant Raman scattering of the OVI resonance lines A 1032 
A and A 1037 A by hydrogen, near the Lyman (3 frequency. The 
OVI photons are produced in the region near the hot component 
and then Raman scattered in the neutral region near the red giant. 
The very existence of this line suggest a temperature of > 100000 
K for the ionizing source (Miirset & Nussbaumer 1994) and a scat- 
tering region with Nh i ^ 10^^ cm~^. The A 6830 A line appears 
as a broad emission in the 2001 spectra of AECir with two peaks 
on the red side at A 6829.2 A and A 6836.6 A (Fig. 7). The peaks 
surround a central depression at A 6831.8 A. Some narrow absorp- 
tions are also observed in the profile, and they could be part of the 
spectrum of the cool star. A single Gaussian fit to the overall A 6830 
A emission line gives EW = 4.66 A, GFWHM = 23.3 A, and Ac 
= 6829.32 A. The emission line central wavelength compares well 
with the mean value of Ac = 6829.5 ± 0.5 A for 60 symbiotic stars 
given by Allen (1980). The shape of the line suggests the presence 
of many components, and the emission wings are extended, as in 
the case of Ha, but reaching velocities up to ~ 1000 km s"'^. 
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3.6 The AA 4640-4650 emission and possible X-ray excitation 

Emission around A 4640 A lias been detected in X-ray sources 
and also in non X-ray emitters like novae, Wolf-Rayet stars, plane- 
tary nebulae, gaseous nebulae, symbiotic stars. Of, Be, and P Cygni 
stars. McClintock, Canizares and Tarter (1975) suggested that this 
emission comes from ions of CIII, NIII and OIL They also dis- 
cussed three possible origins for this emission; non-selective emis- 
sion by one or more of the ions C III, N III and O II, selective emis- 
sion of NIII A A 4634-4641 in the "Of star process" and selective 
emission of NIII A A 4634-4641 in the Bowen fluorescence pro- 
cess. They argue that it appears likely that the Bowen fluorescence 
mechanism is an important source for N III AA 4634-4641 emission 
in X-ray sources. 

McClintock, Canizares and Tarter (1975) list a series of Cin, 
N ni and O II lines that should be detected in the case that non- 
selective emission is the predominant mechanism at work (their 
Table 4). The fact that we do not observe many of the strong lines 
as cm 4186.9, 4325.6, 5695.9 and NIII 4379.1, 4514.9, argues 
against this mechanism. The selective emission of NIII A A 4634- 
4641 in the "Of star process" is also discarded since it should pro- 
duce absorptions at NIII AA 4097-4103, which are not observed. 
Therefore our observations favor the Bowen fluorescence mecha- 
nism for the AA 4640-4650 emission in ae, as usual in symbiotic 
stars. If this view is correct, then we should detect many OIII emis- 
sion lines blueward 3450 A (McClintock, Canizares and Tarter, 
1975). AECir has not been hitherto observed in this region. 

The emission line ratio NIII 4641/4635 has been used as a 
diagnostic for electron density in hot plasmas produced in pulsed 
laser deposition experiments for thin film depositions (Bohigas, 
Perez-Tijerina & Machorro 2004). These authors determined popu- 
lation levels in a 10.000 K plasma considering radiative decay and 
coUisional excitation and de-excitation. Their Fig. 1 indicate that 
in the years 2000 and 2007, when NIII 4641/4635 was between 
0.7-1.1, the electron density was ~ 10^" or ~ 10" cm • (the 
case around 10^ cm~^ is excluded due to the absence of forbid- 
den lines). On the other hand, in 2001 this ratio increased to 2.5 
(Fig. 7 and Table 2). This value is out of the range predicted by the 
aforementioned calculations, but an extrapolation of them implies 
an increase of electron density between both epochs. 

It is possible that the nebula in AE Cir is excited by X-rays, 
the existence of ionized helium argues in favor of this view. AE Cir 
was covered during the ROSAT all-sky survey in August 1990 for 
about 550 sec (Aug 15-17 1990, HJD 2448118.9-2448120.9), but 
no soft X-ray emission was detected. The object was not observed 
by XMM and Chandra, and does not appear in the Galex UV sur- 
vey. The non-detection by ROSAT implies an upper limit of 0.02 
counts s~^ in the PSPC-instrument. Assuming a blackbody spec- 
trum with 50 (200) eV and E{B - V)= 0.3, this implies an unab- 
sorbed flux limit of 3 x 10"" (5 x lO'^^) erg cm s . Assuming 
a thermal bremsstrahlung spectrum with IkeV temperature gives 
7.5 X 10~^^ erg cm~^ s~^ and a power law with photon index 1.0, 
1.5 or 2.0 gives 5, 6, or 8 x 10^ ^'^ erg cm^^ s^^. This corresponds 
to a luminosity upper limit of 3.5 x lO'^'^ - 6 x lO'^'^ (D/10 kpc)^ erg 
s ^ ^ . Note that if the distance were ~ 1 kpc and the X-ray emission 
soft, foreground absorption would make it particularly diffictilt to 
detect X-rays from AE Cir. 



3.7 On the He I emission and Balmer decrements 

The He I emission has been used as a diagnostic tool for the physi- 
cal conditions in the gaseous envelopes of symbiotic stars. We com- 



Table 4. Observed Balmer decrements compared with predictions for 
Case-B recombination theory for different temperatures (Osterbrock & 
Ferland 2006). Typical error for AECir emission line ratios are 10%. 
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0.36 
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0.25 
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2.36 


0.40 


0.37 


0.45 


2500 K 


3.30 


0.44 


0.24 




5.000 K 


3.05 


0.45 


0.25 


0.56 


10.000 K 


2.87 


0.47 


0.26 


0.62 



pared our He I emission line data with models by Almog & Netzer 
(1989) and Proga, Mikolajewska and Kenyon (1994), finding that 
they are in the upper limit of symbiotic He I ratios, in the region of 
S-type symbiotics, and that they are not exactly reproduced by the 
models. However, if we extrapolate the model results to AE Cir, the 
He I emission forming region would has very high optical depth, 
electron density and temperature at the low state. We note also that 
our He I ratios present significant deviations from Case B predic- 
tions, as usual in symbiotic stars. 

We also compared our emission line strengths with predictions 
of the model for an illuminated red giant wind by Proga, Kenyon & 
Raymond (1998). These authors calculated equivalent widths and 
emission line ratios for the emergent spectrum of a red giant wind 
illuminated by the hot component of a symbiotic binary system. 
Different wind velocity laws, mass loss rates and hot component 
temperature were used in the models. A comparison with the 
equivalent widths of AECir shows that the H/3 equivalent width 
is reproduced by a strong illumination model with log Th decreas- 
ing from ~ 4.8 in the low state to ~ 4.6 in the high state, but the 
He I 5875 and Hell 4688 lines require logT^ « 5.2. This could 
be interpreted in principle as two different forming regions for hy- 
drogen and helium lines, but the predicted H/3 equivalent width for 
logTft ~ 5.2 is ^ 100 A, something that is not observed. In ad- 
dition, we found that our He I emission line ratios are larger than 
expected. The failure of these models with the He I A6678/A5875 
emission line ratio is not restricted to AE Cir, but was noted to oc- 
cur with symbiotic stars in general by Proga, Kenyon & Raymond 
(1998). 

Our Balmer decrements are listed in Table 3. The general dis- 
crepancy with Case-B recombination theory probably indicate high 
densities and the importance of atomic colUsional excitation in the 
nebula. When comparing our decrements with those predicted by 
Drake & Ulrich (1980), we found that the 2007 data indicate an 
increase of the electron density at the high state. This qualitative 
result is also obtained when comparing our Balmer decrements 
with those predicted by the models for an irradiated red giant at- 
mosphere by Schwank, Schmutz and Nussbaumer (1997). The ir- 
radiated wind model of Proga, Kenyon & Raymond (1998), on the 
other hand, can reproduce the Ha/H/3 decrement observed in the 
high state only with a model of high irradiation Itmiinosity. 



4 GENERAL DISCUSSION 

As mentioned above, the presence of emission lines, the Raman 
scattered line at A 6830 A and the detection of a cool stellar com- 
ponent indicate a symbiotic classification for AECir. In a sym- 
biotic system, a red giant faces a hot component, normally a hot 
white dwarf, that ionizes part of the red giant wind. UsuaUy symbi- 
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otic star spectra reveal the red giant absorption lines, the emission 
lines coming from a neutral region in the wind and also the high 
ionization species formed closer to the hot component. The rather 
early spectral type of ~ K5, the hot ionizing source and the high 
nebular electron density indicate that AE Cir is a new member of 
the small group of s-type yellow symbiotic stars (Schmid & Nuss- 
baumer 1993). 

The behavior in Fig. 1 suggests the presence of recurrent out- 
bursts in AE Cir. We discarded total eclipses as explanation for the 
low states since traces of both components are observed in the high 
and low states. However, the outbursts are not reminiscent of typ- 
ical symbiotic novae outbursts. While symbiotic novae have out- 
burst amplitudes of a few magnitudes, the rise time is of the order 
of years, and the decay time of order of decades. In the case of 
AE Cir, the rise time in 1989 (HJD ~2447700) could be only 2 
weeks (if we trust in the two visual observations discussed previ- 
ously) and the next rise happened already in 2001, 12 years after 
the previous one. Compared to all other known symbiotic novae, 
AECir would be on the very short-timescale end. The outbursts 
observed in AE Cir can be better classified as of the type called 
"classical symbiotic outburst" (Sokoloski 2003). During a classi- 
cal symbiotic outburst, the optical brightness typically increases by 
several magnitudes, but in some cases can increase by as little as 
one magnitude. The system may take weeks or months to reach 
maximum brightness, and then months or years to fade. However 
a peculiar feature of the AE Cir light curve is the long time spend 
at maximum: the low state has ~ 38% the length of the high state. 
The general cause for classical symbiotic outbursts has not yet been 
established (Sokoloski 2003), although in Z And one outburst was 
identified starting with characteristics of a dwarf nova outburst (i.e. 
triggered by a disc instability) and followed by a nova-type outburst 
(i.e. probably caused by enhanced nuclear burning on the surface of 
the white dwarf, Sokoloski et al. 2006). Incidentally, the shape of 
the main outburst of ae is remarkably similar to a superoutburst of 
a SU UMa type, dwarf nova, which is explained by a thermal-tidal 
disc instability (e.g. Warner 1995). 

We have calculated the temperature of the hot component 
using the relation of Murset & Nussbaumer (1994) that uses the 
ionization potential of the highest-ionization species visible in the 
spectrum Xmax ■ This relation is Th = 1000 Xmax and gives 100 000 
K in the low state in 2000 and 2001 and 55 000 K in the high state 
in 2007. We also used the lijima (1981) relationship for a radiation- 
bounded nebula, based on the ratio of nebular emission Unes. We 
obtained upper limits for Th of 208 000 K, 268 000 K and 242 000 
K in 2000, 2001 and 2007. The lower calculated temperature of the 
hot component in the high state is deduced from the disappearance 
of the 6830 A emission feature. We observe that the decrease of Th 
argues against the accretion disc instability model for the outbursts 
of AE Cir. In this model the disc increases luminosity during out- 
burst due to an increase of mass accretion rate (e.g. Duschl 1986), 
but in this case we should observe an increase of the hot compo- 
nent temperature. On the other hand, the nuclear burning model 
imphes an expansion of the white dwarf photosphere at constant 
bolometric luminosity and predicts a decrease of the hot compo- 
nent temperature, as observed. However, the very fast development 
of the main outburst is contrary to the model predictions; the shell 
envelope should expand on the thermal time scale, which is years to 
decades for typical conditions (Iben 1982, Iben & Tutukov 1996). 
The second outburst instead has a much more gradual development, 
more in agreement with the nuclear burning model. However, we do 
not observe spectral features that could reveal an expanded white 
dwarf atmosphere, only the presence of a blue continuum at the 



high state (Fig. 5). In particular, no A-type or F-type shell spec- 
trum is observed as in the symbiotic stars CI Cyg and RR Tel at the 
high state (Belyakina 1979, Thackeray 1950). A K-band spectrum 
of AE Cir kindly taken by Elisa Nespoli (private communication) 
with the ESO NTT SOFI in May 25, 2007 (MJD= 2454245.9942) 
with a S/N of 35 shows no Br7 emission or absorption, indicating 
that the nebular Unes are not prominent in the infrared at the high 
state, and that there are no signs of an hypothetical hot stellar com- 
ponent. The red giant is not detected in this spectrum, possibly due 
to the low S/N of the data. 

If full visibility of the red giant occurs at the low state at V ~ 
15.5, this would indicate a distance of 9.4 kpc, for a typical K5 
giant behind the galactic disc. This should place s below 1.4 kpc 
from the galactic plane. We also note that if the high state is due 
to the hot component outburst, its V ~ 13 mag combined with 
log Th > 4.6 should imply a very high Mboi > -6.4, above the 
Eddington limit for a t5T)ical 0.6 M© symbiotic white dwarf, but 
around the Eddington limit for a massive white dwarf accretor of 
~ 1 Mq (Belczynski et al. 2000, Mikolajewska 2003). 

In summary, while the behavior of Th and the long outburst 
rising time could suggest nuclear burning, the blue continuum and 
the lack of A-type or F-type shell spectrum at the high state point 
to the accretion disc instability model. We now examine if the ex- 
planation for the outbursts of AE Cir could be related to the bizarre 
behavior of the spectral features of the secondary star. 

It is remarkable that we observe the TiO band around A 6300 
A at roughly comparable depth in all spectra of Fig. 5. The differ- 
ence in V magnitude between the 2007 and 2000 spectra is about 
1.5 magnitudes. This should imply a change of the relative contri- 
bution of the cool giant to the total light by a factor 4, but the change 
in equivalent width of the TiO band is only 25%. In other words, 
it should be practically impossible to detect the secondary star at 
maximum if its brightness remains the same that at minimum. The 
fact that we observe the secondary star contributing roughly the 
same fraction of light during the outburst cycle suggests that: (i) 
the star varies its brightness or (ii) the whole system light is ob- 
scured by some source during the low state. At the same time, the 
spectral type determinations indicate that color and temperature of 
the giant star remained roughly unchanged during the observation 
epochs. 

From Fig. 1 the overall brightness at the epoch of the 
SMARTS observation is V 12.7 mag. If we assume that half of 
the luminosity comes from the giant star we get V ~ 13.5 for the 
giant star. Using this value and E{B — V)= 0.26 we obtain a dis- 
tance of 10.2 kpc, in nice agreement with our earlier estimate. From 
the above it seems that the putative change of the stellar brightness 
of the giant star in AECir could be large. A ~ 2 mag change at 
constant radius should imply a change in temperature by a factor of 
1 .6, which is not observed. This means that the giant star in AE Cir 
should change its radius by a factor ~ 2, likely engulfing their com- 
panion. Such a big change in radius is hard to understand. 

Regarding the alternative of obscuration, it is compatible with 
the larger reddening observed at minimum, as derived from the 
strength of the NaD doublet and system color, but not with the 
larger emission observed at minimum neither with the long de- 
crease of brightness observed after maximum. We conclude that the 
data present difficulties for the canonical symbiotic outburst mod- 
els. 

We note that if the ellipsoidal like variability with a pe- 
riodicity of 342 days is the binary period, the amplitude, if 
ellipsoidal in nature, largely exceeds the expected values for a 
tidally distorted star (Soszynski et al. 2004, Morris 1985). We note 
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that the B — R color shows a systematic change over the 342 
day period (see Fig. 3), being maximum at the time of minimum 
flux of the 342 day oscillation. This would be expected if the 
cool star is in front of the hot component. If this were true, then 
the rather narrow eclipse light curve observed in Fig. 3 could 
indicate that the red giant self-eclipses a zone of chromospheric 
fluorescence, as proposed for VI 329 Cyg (Chochol & Wilson 
2001). The visibility of the cool giant at minimimi and max- 
imum could in principle allow detection of such a variability. 
If we really observe a chromospheric eclipse in the SMARTS 
data, the 40-day eclipse length suggests a radius for the red giant of: 

^9 = 7 X si X 27ra = °-f^a (2) 

where a is the binary separation and / a geometric factor depend- 
ing on the characteristics of the illuminated surface. If / = 2, i.e. 
the eclipse corresponding length is comparable to the stellar diam- 
eter, the radius is surprisingly close to the Roche-lobe radius for the 
secondary star in a binary of mass ratio unity, viz. 0.38a (Warner 
1995). As the mass of a K5 giant is ~ 1.2 Mq and the mass of a 
white dwarf (assumed to be the hot component) is around 1 Mq, 
then it is possible that the red giant in AECir fills its Roche-lobe. 
If we use the above masses for the components and the third Ke- 
pler law we obtain a = 1.87 x 10^^ cm and Rg = 6.92 x 10^^ 
cm = 99 R0 . The condition of Roche-lobe overfilling, along with 
the putative orbital period of 342 days, constrains the stellar mean 
density to logp = -5.8 g cm~^ (e.g. Warner 1995). This value yield 
an independent estimate for the red giant mass of 1.1 Mq, in nice 
agreement with our early assumption. We note that the derived ra- 
dius is in the upper limit of radius for red giants with a temperature 
of 3950 K (van Belle et al. 1999). Models for irradiated late type 
giants show an expansion of the atmosphere by several percents 
for large hot component luminosities (Proga et al. 1996), hence ir- 
radiation could play a role in the size of the red giant in AE Cir. 
Irradiation is also suggested by the emission line ratios discussed 
in the previous section. The above calculation implies a luminos- 
ity for the red star of Lboi = 7.80 x 10^** erg s"\ yielding Mboi 
= -3.42 and using a bolometric correction of BC = -1.02 for a K5 
giant suggests My = -2.4. 



5 CONCLUSIONS 

In this paper we have investigated the nature of the symbiotic star 
candidate AE Cir. We have analyzed new optical photometric and 
spectroscopic data, 2MASS infrared photometry and 24 years of 
visual photometry in an integrated way in order to get a better 
understanding of the system. Our conclusions can be summarized 
as follows: 

• The symbiotic nature is confirmed. This result is based on the 
detection of the A 6830 A emission feature and the spectral signa- 
tures of a cool stellar component of spectral type ~ K5. 

• The spectral type ~ K5 and the spectral and photometric prop- 
erties indicate that ae is a member of the small group of yellow 
symbiotic stars of the s-type. 

• The fight curve is characterized by outbursts lasting ~ 4000 
days and overall amplitude of variability about 4 magnitudes. The 
outbursts show a slow decline of ~ 2 x 10"'' mag/day. The duration 
of the low state is about 38% the high state. 

• A strong signal at 342 ± 15 days is detected in the fight curve. 



The light curve folded with this period shows two broad minima 
with different amplitude. 

• The spectrum in the low state is characteristic of a symbiotic 
star without forbidden lines and very strong He 1 4868 emission. At 
maximum the emission feature at A 6380 A disappears, the spec- 
trum becomes flat and the relative intensity of the He 1 4686 A line, 
compared with the Bafiner lines, becomes lower. At the same time 
the overall emission fine spectrum shows smaUer equivalent widths. 
Signatures of the cool companion are observed at the high and low 
state. 

• We observe H I, He I and He II emission line variability in time 
scales of minutes in 2001 being larger in the line centers. The NaD 
and A 6820 A lines do not follow this variability. We also observe 
larger broadening and asymmetry in the Balmer lines compared 
with others emission fines in this epoch. 

• The data suggest that full visibility of the red giant should 
occur at the low state at V ~ 15.5. For a typical K5 giant behind 
the galactic disc this would indicate a distance of 9.4 kpc. This 
should indicate that se is about 1.4 kpc below the galactic plane. 

• We suggest that an eclipse-like event could be interpreted as 
self-occultation of the irradiated red giant hemisphere. In this case 
we found that it is possible that the red giant fiUs its Roche-lobe. 
Assuming a white dwarf hot component of 1 Mq, we estimate a 
radius of 99 Rq , a mass of 1 . 1 Mq and Mv = -2.4 for the red giant 
at the epoch of the SMARTS observations. 

• The atypical character of se is revealed in the long time pass- 
ing in outburst and the dwarf nova superoutburst shape of the light 
curve. Another atypical feature is the rather strong TiO absorption 
band around A 6300 A at the high state. This could indicate that the 
red giant and the hot component are partially obscured at minimum 
or the brightness of the secondary star is not constant. Obscuration 
hardly explains the large line emission observed at minimum and 
the long brightness decay after maximum whereas variations of the 
red giant brightness imply unrealistic changes in the stellar radius. 

• At present, the available data for AE Cir are hard to reconcili- 
ate with canonical models for symbiotic outbursts. 
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